Exact ground state and elementary excitations of the spin tetrahedron chain 
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We study the antiferromagnetic spin exchange models with S = 1/2 and S' = 1 on a one- 
dimensional tetrahedron chain by both analytical and numerical approaches. The system is shown 
to be effectively mapped to a decoupled spin chain in the regime of strong rung coupling, and a spin 
sawtooth lattice in the regime of weak rung coupling with spin 2S on the top row and spin S on 
the lower row. The ground state for the homogeneous tetrahedron chain is found to fall into in the 
regime of strong rung coupling. As a result, the elementary excitation for the spin- 1/2 system is 
gapless whereas the excitation for the spin-1 system has a finite spin gap. With the aid of the exact 
diagonalization method, we determine the phase diagram numerically and find the existence of an 
additional phase in the intermediate regime. This phase is doubly degenerate and is characterized 
by an alternating distribution of rung singlet and rung spin 25. We also show that the SU{<i) 
exchange model on the same lattice has completely different kind of ground state from that of its 
SU{2) correspondence and calculate its ground state and elementary excitation analytically. 

PACS numbers: 75.10.Jm 
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I. INTRODUCTION 



The study of quantum antiferromagnetic spin mod- 
els with strong frustration has attracted great atten- 
tions over the past decades. Early investigations of the 
frustrated quantum magnets were partly motivated by 
the work of Anderson to search the resonating-valence- 
bond (RVB) ground state in such systems^. In the frus- 
trated quantum magnets, the magnetic ordering is gen- 
erally suppressed by the frustration. Some well-studied 
frustrated magnetic systems include, for example, the 
kagome lattice and pyrochlore lattice, in which the inter- 
play between frustration and quantum fluctuation leads 
to rich varieties of phenomena. Recently, the frustrated 
magnets are believed to be prominent candidates of re- 
alizing spin liquid states with exotic ground state and 
deconfined fractional excitations^''^'"'. While the mecha- 
nism of deconfinment in two-dimensional (2D) magnetic 
systems is less clear, the deconfinement of spinous in 
quasi-one-dimensional magnetic system, which is closely 
related to the phenomenon of spin charge separation, 
is well investigated^iii^. Since 1980's, a large number 
of low-dimensional frustrated magnets have been syn- 
thesized experimentally-. Generally speaking, strongly 
geometrical frustration in these systems allows the sim- 
ple dimerized state to be the ground state of the low- 
dimensional frustrated spin system and opens a spin gap. 
So far, a variety of quasi-one-dimensional frustrated mod- 
els have been studied theoreticall y^°i^^'^^ . Additionally, 
important progress has been made in trapping cold atoms 
under a highly controllable way very recently, and thus 
it stimulates intensive investigation on how to simulate 
the magnetic systems using cold atoms. A number of 



schemes have been proposed to implement a variety of 
quantum spin models in optical lattice&iS. 

In this article, we investigate both the ground- and 
excited-state properties of the spin models on a one- 
dimensional (ID) tetrahedron chain as shown in Fig. 1, 
with site spins Si residing in four corners of each tetra- 
hedron. The basic unit, i.e, a tetrahedron is composed 
of four spins with equal antiferromagnetic exchanges be- 
tween each pair of spins. The tetrahedron chain can be 
also viewed as a ID pyrochlore strip in a 2D pyrochlore 
lattice, in which only two of four corners of each tetrahe- 
dron are shared by a neighboring tetrahedron. Generally, 
a three-dimensional (3D) pyrochlore lattice is a network 
of corner-sharing tetrahedra and a 2D pyrochlore model, 
named also as a checkerboard-lattice model, is obtained 
by a projection of the 3D lattice on a plane. Different 
from the ID pyrochlore strip considering in this article, 
for both the 3D and 2D pyrochlore lattices, each corner 
of the tetrahedron is shared by a neighboring tetrahe- 
dron. As one of the most frustrated antiferromagnets, the 
model of spin pyrochlore lattice has been investigated by 
a variety of techniques including the semiclassical large-5' 
limit, large- iV expansion of the SU{N) model, the con- 
tractor renormalization method based on the cluster ex- 
pansion, and the bosonization method on the anisotropic 
limii*i^ii^ii^iii. In spite of the intensive research, even 
the ground state properties of the 3D pyrochlore lattice 
are not well understood. For the 2D pyrochlore lattice, 
the numerical results based on exact diagonizations have 
shown that the ground state has plaquette order—. How- 
ever, for the ID pyrochlore strip, we can determine its 
ground state and elementary excitation in an exact man- 
ner. With the aid of numerical diagonalization of the cor- 
responding spin lattice systems with small sizes, we also 
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FIG. 1: (color online) Schematic pictures of (a) the homoge- 
neous spin tetrahedron chain, (b) the spin tetrahedron chain 
corresponding to model (|2]). 



investigate the quantum phase transitions of the ground 
state due to the change of exchange strengths along per- 
pendicular rungs. 

The spin model on a ID pyrochlore strip as shown in 
Fig.l is described by the Hamiltonian: 



H — J ^ Si ■ Sj, 



(1) 



where (ij) denotes sum over all the nearest neighbors 
along the tetrahedron chain and Si represents the spin 
operator residing in site i. In this work, we study both 
the spin-1/2 and spin-1 models on the pyrochlore strip. 
Our results show that the excitation spectrum for the 
spin-1/2 system is gapless and the elementary excitation 
of the spin-1 system has a finite spin gap. This model can 
be extended to the cases where the strengths of bonds 
among the tetrahedra are not homogeneous. Here we 
only consider the inhomogeneous case as shown in Fig. 
lb, where we use J± to represent the exchange strength 
along the vertical rung of each tetrahedron. For conve- 
nience, we rewrite the Hamiltonian of the spin tetrahe- 
dron chain corresponding to Fig. lb as 

Ha = j'^2^Si^i{Si^2 + Si^s) + {Si^i + Si^2 + Si^3)Si+i^i 

i 



It is obvious that the model ^ reduces to the model 
(1) when J± = J, i.e., the homogeneous tetrahedron 
chain (1) is a special case of the model Since the 
model can be represented as a sum of local Hamiltonian 
on each tetrahedron, a classical ground state is obtained 
whenever the total spin in the tetrahedron is zero for a 
homogeneous model. It is straightforward that the clas- 
sical ground states have a continuous local degeneracy. 



II. SPIN-1/2 LATTICE 

Firstly, we consider the spin 1/2 case. As we have 
mentioned in the introduction, there is a fundamental 
different property between the ID pyrochlore strip and 
its high-dimensional analogies: not all spins on a tetra- 
hedron are equivalent. The spins on the vertical rung of 
the tetrahedron are not coupled to the neighboring tetra- 
hedron. We define the total spin on each vertical rung as 
Ti = Si^2 + Si^3. It is obvious that is conserved. This 
property enables us to simplify the ID model greatly. 

For the tetrahedron chain model (2), it is convenient 
to reformulate Hamiltonian as 



Ha — J Si^iSj-i-i^i + J (Sj^i + Sj-i 



■T, 



J2T^~NJ^SiS+l)., 



(3) 



where = (T; + 1) with T; = 0, • • • , 2S'. In the strong 
coupling limit Jj^ — > cx), a pair of spins on each rung 
would form a singlet (spin dimer) with Ti = 0. This 
implies that the spins in the horizontal direction along 
the chain is effectively decoupled with the spins on the 
vertical rungs. Therefore, the ground state of Ha is a 
product of the ground state of spin chain and rung sin- 
glets. Explicitly, it is represented as 



\GS) = \BA)^Y[[S,^2,S,,3] 



(4) 



where \BA) denotes the Bethe-ansatz ground 
state wavefunctions of the ID Heisenberg chain 

Hchain = J ^i Si^iSi+i,ia,nd [Si^2,Si,3] = 

([t].,2 [i]»,3 - [i]^,2 [T].,3) /V2 is the dimer singlet 
across the ith vertical rung. The corresponding ground 
state energy is 



(N) - -NJ^, 



(5) 



where {N) is the ground state energy of the A^-site 
Heisenberg spin-1/2 chain. From the Bethe-ansatz solu- 
tion of Heisenberg model, we know the exact ground state 
energy E^^ (N) /N = -0.4431J at the infinite length 
limit. 

In fact, utilizing the Raleigh- Ritz variational 
principleiiii2i'^i2i, we can exactly prove the state 
given by Eq. ^ is the ground state of Hamiltonian ^ 
as long as J± > 2 J. To see it clearly, we can rewrite 
the Hamiltonian (2) with J± = 2J as the sum of a 
Heisenberg chain and 2N projection operators, which 
reads 

Ha — J 5^45*^4.1,1 + —P^^'^{Si^i, Si^2,Si^3) 



P^^^('^«^2, 5*1,3, 5*^+1,1) — 1 



2 L 



(6) 



3 



where 



P^^'^ {Si^i, Si^2, Si^s) — - 



(7) 



is a projection operator which projects a three-spin state 
composed of S'i.i, Si^2 and S^^a into a subspace with to- 
tal spin 3/2. Now it is clear that the state given by Eq. 
(2]) is the ground state of the global Hamiltonian because 
it is simultaneously the ground state of each local sub- 
Hamiltonianiii22iiSi. With the same reasoning, the state 
^ is of course the ground state of Hamiltonian Eq. (2) 
for J± > 2J. It is not hard to check that the state ^ 
is an eigenstatc of the global Hamiltonian Ha by uti- 
lizing the identities (^Si^2 + S'i.s^ [S'i,2, 'S'i,3] = 0, however 
such an eigenstate is not necessary the ground state for 
arbitrary J±. We note that the condition J± > 2J for 
the existence of the ground state given by ) is just a 
sufficient condition which is a very strong restriction. In 
fact, it can be released to a wider parameter regime. We 
expect that there is a critical value and the system 
evolves into another quantum ground state when J± is 
smaller than J^. 

For the homogeneous point J± = J which we are par- 
ticularly interested in, although the above proof is no 
longer applicable, we can still argue that Eq. (|4|) remains 
the ground state, and prove this result by using the nu- 
meric exact diagonalization method. With the aid of 
numerical diagonalization, we may determine the phase 
boundary of the model ([2]) precisely. Since the total spin 
in every rung is conserved, the vertical rungs are either in 
singlets or triplets for a spin-1/2 pyrochlore strip. There- 
fore, the relevant eigenstate of the pyrochlore strip can 
be classified by the values of the total spins on the ver- 



tical rung! 



.22.23 



For convenience, we use E {N, M) to 



represent the eigenenergy of the state with N — M spin 
singlets and M triplets on the vertical rungs. It follows 
that the eigenenergy is given by 



E (iV, M) = -^NJx_ + MJx_ + JE^/2,1 {N, M) , (8) 

where £'1/2,1 {N,M) represents the energy of the lattice 
composed of N spins with S = 1/2 on the sites of the 
lower row and M spins with T = 1 on the top row. 
For each class of state with N — M spin singlets and M 
triplets on the vertical rungs, there are altogether 
different configurations. 

It is instructive to firstly consider the following two 
kinds of configuration which correspond to two opposite 
limits of the vertical exchange: (1) all the states on the 
vertical rungs are spin singlets and (2) all the states on 
the vertical rungs are spin triplets. In the first case, the 
original model can be mapped to a spin chain model as 
displayed in Fig. 2a and the corresponding eigenenergy 
is given by 







(b) 




E{N,Q) 



-NJ^ + JEy2 {N) , 



(9) 



(C) 



FIG. 2: Three types of ground state for the spin-1/2 and 
spin-1 pyrochlore strips. The open and close circles denote 
the spins with S = 1 and 5=1/2 respectively for the spin- 
1/2 pyrochlore strip. For spin-1 pyrochlore strip, the open 
circles represent the spins with S = 2 and the close circles 
represent the spins with S = 1. 



where £1/2 {N) = -Ei/2,1 {N, 0) is the energy of the A^-site 
Heisenberg spin-1/2 chain and the Bethe-ansatz ground 
state energy E^^^ (N) /N = -0.4431 as iV ^ 00 . In 
the second case, the tetrahedron chain can be effectively 
described by a A-chain model consisted of N spins with 
5 = 1/2 on the sites of the lower row and N spins with 
T = 1 on the sites of top row as shown in Fig. 2c and 
the eigenenergy can be represented as 

E{N,N)^^NJ,_ + JEi/2,iiN,N), (10) 

where £1/2,1 {N,N) denotes the eigenenergy of the cor- 
responding A-chain model. The mixed spin A chain can 
be viewed as an alternating spin-l/2-spin-l chain with 
an additional next-nearest-neighbor interaction between 
the spins with Si = 1/2. It is well known that there 
is a ferrimagnetic long-range order in the ground state 
of quantum ferrimagnetic Heisenberg chain^iS^i^. The 
additional interaction between the spins with S' = 1/2 
is a kind of frustration which makes it harder to com- 
pensate the spin with S = 1, therefore the long-range 
order still exists. Although no exact analytical results 
for the mixed spin A chain are known, its ground state 
energy may be determined by numerical exact diago- 
nalization method. For a mixed spin A chain with 
size of 8 -|- 8, we get its ground state energy given by 
£f/2 1 (8' 8) = -0.646773. We note that the ground 
state of a spin A chain or a spin sawtooth model with 
S^ = l/2 is exactly knowni^i^i^i. 

It is clear that the state with all singlets or triplets on 
the rungs are the ground state of the tetrahedron chain 
in the limit of J_l — *■ 00 and J± —^ —00 respectively. 
Flipping a rung singlet into triplet costs an energy of 
Jl , therefore the effect of the antiferromagnetic coupling 
Jl is to prevent the spins on the rung forming triplet. 



4 



Ru^SiHS-tifh-lfidiii 
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FIG. 3: The phase diagram for spin-1/2 tetrahedron chain 
with variable vertical exchange J± in the parameter space of 
c = J_l/J- For c > ci, the system is in a decoupled phase 
whose ground state is a product of rung singlets and the crit- 
ical spin liquid phase on the horizontal spin-1/2 chain; For 
C2 < c < Ci, the ground state is a double degenerate mixed 
state with the alternating spin singlet and spin triplet on the 
rungs; For c < C2, the ground state is a ferrimagnetic state 
with spin triplet on the rungs. 



On the other hand, a triplet in the rung acts effectively 
as a spin with 5* = 1 which lowers the total energy by 
interacting with its neighboring spins with S — 1/2 on 
the lower rung. The competition between the two pro- 
cesses gives rise to the complexity of the phase diagram 
for the tetrahedron chain. An interesting question aris- 
ing here is whether some intermediate phases exist be- 
tween the phases with fully paired singlets and triplets 
on the rungs. To determine the phase boundary of ^ nu- 
merically, in principle, we need consider all the different 
configurations of Tj on the rungs. Among a given class, 
we find that the configuration with the spins on the top 
row repelling each other has lower energy. For example, 
as shown in Fig. 2b, the configuration with alternating 
spin 1 and spin on the top row has the lowest energy 
among the Cj^^^ configurations. After considering all the 
rung configurations, we get a phase diagram as shown in 
the Fig. 3. As expected, there is an intermediate phase 
which is effectively described by the ground state of its 
equivalent model as shown in Fig. 2b. Corresponding to 
Fig. 2b, there is another equivalent configuration which 
is obtained by totally shifting the spins on the top row a 
lattice space. In this phase, the triplet and singlet on the 
rungs distribute in an alternating way and the ground 
state is doubly degenerate. 

Next we will give a description on how to determine 
the phase diagram as shown in Fig. 3. Since the origi- 
nal pyrochlore strip can be classified by the values of the 
total spins on the vertical rungs or equivalently by M , la- 
belling the numbers of the rung triplets, the ground-state 
energy of our tetrahedral spin chain, for a given value of J 
and Jj_ , will be given by the lowest value of Eq. ([8]) with 
A/ = 0, 1, • • • , A^. In Fig. 4, we show as a function of 
c = J_l/J all the lowest eigenenergies for the lattice sizes 
with 4, 6, and 8 tetrahedra. The energies corresponding 
to the point of c = (intercept points with the verti- 
cal axis) are the lowest energies of £'1/2,1 (A^, -Af), from 
above to below, with M = 0, 1, • • • , A^. The slope values 
of the straight lines depend upon the number M . Our 
numerical analysis shows that the ground state in the 
whole parameter space c is determined by three kinds 
of configurations with M = 0, A^/2, A^, corresponding to 
the configurations of (a), (b) and (c) displayed in Fig. 2. 
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FIG. 4; (color online) The lowest energies of E{N, M)/J as a 
function of c = J±/ J for = 4, 6, 8. The solid, dashed and 
dotted lines corresponding to the case with M — N , M = N/2 
and M = or (c), (b), (a) as displayed in Fig. 2. 



and their crossing points determine the phase transition 
points. For simplicity, we omit all the excited energy lev- 
els for the case of A^ = 8 in the Fig. 4. From our numer- 
ical results, we obtain the up critical value ci = 0.9529 
and the down critical value C2 = 0.7214 for the original 
pyrochlore strip with a size of 24 sites (8 tetrahedra or 
A^ = 8). Similarly, we get a = 0.88186 and C2 = 0.70546 
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FIG. 5: (color online) The phase boundaries ci and C2 versus 
the sizes of the system. 
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FIG. 6: (color online) The ground state energies of systems 
(a), (b) and (c) versus the sizes of the system. 



for TV = 4, ci = 0.92952 and ci = 0.72351 for iV = 6. 
Here C\ corresponds to the crossing points of lines (a) 
and (b) , whereas C2 corresponds to the crossing points of 
lines (b) and (c). Since we need much smaller memory 
size to diagonalize systems of (a) and (b) than system of 
(c), therefore we can calculate even larger system to de- 
termine ci. For example, we get ci — 0.96228 for iV = 10 
and ci — 0.96752 for iV = 12. In Fig. 5, we analysis the 
finite size scaling of ci and ci. The linear fit of c\ and ci 
gives ci = 0.9767 ±0.0020 and C2 = 0.73021 ±0.00577 for 
N oo. We can also determine the phase boundaries 
ci and C2 in an alternative way. By extrapolating the 
ground state energy of (a) and (b) to infinite size, we then 
determine ci and C2 by the crossing of the energy levels 
of (a), (b) and (b), (c) respectively. In Fig. 6, we display 
the finite size scaling of the ground state energies of (a), 
(b) and (c). The ground state energies per site, obtained 
by linear fit to the infinite limit, corresponding to the 
configurations of (a), (b) and (c) are 0.44261 ± 0.00019, 
0.62033 ±0.00091 and 0.64675 ±0.00001 respectively. By 
this way and using Eq. ([5]), we get ci — 0.9757 and 
C2 = 0.72601 in the limit of infinite size. If we use the 
ground state energy obtained by Bethe ansatz for (a), we 
get ci = 0.9748. 

Before ending the discussion of the spin- 1/2 pyrochlore 
strip, we would like to give some remarks to a general- 
ized spin pyrochlore strip where the horizontal exchange 
is variable. With the same reasoning as that of the model 
([2) , we can easily get the sufhcient condition for the exis- 
tence of the fully dimerized state on all the vertical rungs, 
which reads J± > 2 J and is irrelevant with the horizontal 
exchange. When the horizontal exchanges are zero, the 
model has totally different ground sates and it falls into 
the class of ID diamond mod o^^i^^ , which has also fully 
dimerized state on the vertical rungs as the ground state 
but the ground state is highly degenerate with a degen- 
eracy of 2^ because the unpaired N spins are completely 
free. 



III. SPIN-1 LATTICE WITH SU(2) SYMMETRY 

The above method can be directly extended to deal 
with the spin-1 case. For the spin tetrahedron chain ^ 
with S = 1, we can prove that the ground state of ([2|) 
is a product of fully dimerized singlets on the rungs and 
the ground state of horizontal spin-1 chain as long as 
J± > 4J. Explicitly,the ground state can be represented 
as 

\GS) = \Haldane) (g) [] [S,,2, ^,,3] (11) 

i 

where \Haldane) denotes the ground state of the ID spin- 
1 chain and 

m — — 1 

is a spin singlet across the ith vertical rung. The proof 
is rather similar to its spin-1/2 correspondence and is 
straightforward when we rewrite the Hamiltonian ([2]) in 
the following form 

Ha — E/ ['^^('^''l' "^i.s) + hA(<S'i,2, 'S'i,3, S'i+1,1) 
i 

+ JY,S^sS^+l,l (12) 

i 

with hA(S'i4, 5'j,2, 5'j,3) = JSi^i{Si^2 + S'i.s) + '^<5'i,2 • 
Si^3 and hA(5'i^27 "Si^a, 5^+1^1) = J{Si^2 + 5*2, 3)5*4+1,1 + 
^Si_2 ■ Si^3. Now it is easy to find that the ground state 
of hA is a product of paired singlet on the rung and 
a free spin on the unpaired site as long as Jj_ ^ 4J. 
By the variational principle, we conclude that the state 
(fTTj) is the ground state of spin-1 model given by Eq. 
(2) for J± ^ 4J. Certainly, the sufficient condition for 
the existence of fully dimerized ground state on the rung 
can be released to a lower bound. In principle, we can 
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FIG. 7: The phase diagram for spin-1 tetrahedron chain. For 
c > Cl, the system is in a decoupled phase whose ground 
state is a product of rung singlets and Haldane phase on the 
horizontal chain; For C2 < c < ci , the ground state is a double 
degenerate mixed state with the alternating spin singlet and 
spin quintet on the rungs; For c < C2, the ground state is a 
ferrimagnetic state with spin quintet on the rungs. 

determine it numerically following a similar scheme of 
the spin- 1/2 case. 

Since the total spin Ti on each vertical rung can be 
0, 1 or 2, the spin-1 tetrahedron chain can be mapped 
to a mixed sawtooth-like model according to its configu- 
rations of the vertical rungs. Several relevant configura- 
tions are: (1) Spin singlets on all the vertical rungs with 
the eigenenergy given by 

Et=q = -27V Jj_ + JEi (TV) , (13) 

where Ei (N) denotes the eigenenergy of a spin-1 chain 
of N sites. From the known numerical results'^*', we get 
ground state energy per site Ef (N) /N = —1.4051 in the 
large N limit. (2) Spin triplets on all the vertical rungs 
with the eigenenergy given by 

Et=i = -NJj_ + JEi^i {N, N) . (14) 

Here Ei^i {N, N) represents the eigenenergy of a spin-1 
A chain of 2N sites^i. (3) Spin quintet (T = 2) on all 
the vertical rungs with the eigenenergy given by 

Et=2 - NJ^_ + JEi^2 {N, N) . (15) 

Here i?i,2 {N, N) represents the eigenenergy of a mixed 
spin-1 and spin-2 A chain with spin 1 on the site of the 
lower row and spin 2 on the site of top row. Just simi- 
lar to the spin-1/2 pyrochlore strip, we need consider all 
the different configurations of the rungs. Numerically, 
we find the phase diagram is rather similar to that of the 
spin-1/2 case. As displayed in Fig. [71 there is also an 
intermediate phase between the fully paired singlet state 
with Ti — and state with all Ti — 2. The intermedi- 
ate state is twofold degenerate with the configuration of 
alternating singlet and quintet on rungs. Such a state 
can be schematized in terms of Fig. 2b with the open 
circle denoting the state with Ti = 2. From our numer- 
ical results, we obtain ci = 0.91596 and C2 = 0.88404 
for the original tetrahedron chain with a size of 18 sites 
(N=6) as well as a = 0.9134 and C2 = 0.8452 for iV = 4. 
With similar reasoning as the spin-1/2 case, we can get 
Cl = 0.94189 for even larger system with N = 8. Linear 
fitting of datum of Ci versus 1/iV^ and extrapolating it to 
the Hmit of the infinite size, we get ci = 0.95983±0. 02748 
for — > oo. 



As a natural generalization, it is straightforward to ex- 
tend the spin pyrochlore model with SU (2) symmetry to 
the case with arbitrary spin S. For the spin-5 model, 
a sufficient condition for the existence of rung-dimerized 
ground state is Jj, > 4 — 3/(s -I- 1) for half-integer spin 
and J_L > 4 for integer spin. In the dimerized phase, the 
horizontal chain is decoupled with the spins on the rungs 
and therefore the elementary excitation of the tetrahe- 
dron chain is gapless for half-integer-spin model or opens 
a gap for the integer-spin model. 



IV. SPIN-1 LATTICE WITH SU(3) SYMMETRY 

For the spin-1 system, the most general Hamiltonian 
has a biquadratic exchange term besides the bilinear term 
and it exhibits much richer quantum phase structures 
than the bilinear model^^. When the biquadratic ex- 
change has the same strength of the bilinear exchange, 
the Hamiltonian H = J^i^ij) with 

hiS,,S,)^S,-S, + {S,-Sjf (16) 

owns the SU{3) symmetry. For spin-1 systems of tran- 
sition metal compounds where two electrons are cou- 
pled ferromagnetically by Hund's rule, the biquadratic 
exchange term originates from a fourth order perturba- 
tion process. Its magnitude is thus small compared to 
the bilinear terms, and thus the SU{3) symmetry is not 
applicable. However, in the cold atomic physics, most 
atoms have high hyperfine multiplets, thus it is possi- 
ble to achieve high symmetries. For example, the ®Li 
atom is with nuclear spin 1 and electron spin 1/2. In 
a weak magnetic field, electron spin is polarized, while 
nuclear spin remains free. Recent studies indicated that 
the three nuclear spin components can be described by 
an approximate SU (3) symmetry^. 

It is well known that the SU{3) exchange model on 
a chain has quite different properties from that of the 
SU{2) bilinear model. Therefore we may expect that the 
SU{3) tetrahedron chain 

3 

H = J±^h{Si^2,Si^3) + J^^h{Si^a,S^+l^l) + 

i i a— 1 

[/i(5,,i, 5,,2) + h{S,,i, ^,,3)] (17) 

i 

also displays different phase structure from that of its 
SU{2) correspondence ^ . The model with J± — J 
was initially proposed by three of us with S. C. Zhang 
in Ref.-- as an example of the SU{N) generalization of 
the Majumdar-Gosh model, however no analytical results 
have been given there. Observing that the Hamiltonian 
can be written as a sum of the Casimir of the total spin in 
each tetrahedron and the representations with the small- 
est Casimir made out of four sites in the fundamental 
representations is three dimensional, we concluded that 
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the state of trimer products is the GS of the SU (3) tetra- 
hedron chain. The ground state is twofold degenerate. 
ExpHcitly, the ground state of the SU{3) spin tetrahe- 
dron chain can be represented as 

\GS),=l[T{S,,uS^,2,S,,^) (18) 

i 

or 

\GS)^ = Y[TiS,,2,S,,3,S,+i.,i), (19) 

i 

where 

T{Si,Sj,Sk) = ;y|'^"/37 \(^i^P3,lk) 

represents a trimer state which is a singlet composed of 
three spins on site i, j, and k . Here ai denotes the spin 
on site i with the value a taking 1, 0, or —1 and ea/s-y is 
an antisymmetric tensor. 

In the following, we shall calculate the ground state 
energy and elementary excitation of Eq. (jl7l ) analyt- 
ically. For convenience, we make a shift of constant 
to the Hamiltonian ^Tf] ) by replacing h{Si,Sj) with 
Pi,j = h{Si,Sj) — 1. We note that modification of J± 
does not lift the degeneracy of the left- and right-trimer 
states. For J± > J, the state of trimer products is of 
course the ground state of the SU(3) tetrahedron chain 
and the corresponding ground state energy is 

Eg = -2NJ- NJ±. 

Breaking a singlet of trimer will cost a finite energy, thus 
the elementary excitation of the SU(3) tetrahedron chain 
has an energy gap. For a three-site cluster, the trimer 
singlet is represented by a Young tableaux [l'^] and the 
first excited state above the singlet are represented by the 
Young tableaux [2^1]. When a trimer singlet is broken, 
it decomposed into a monomer and a paired dimer. For a 
system with degenerate ground state, the monomer and 
dimer can propagate freely in the background of trimer- 
ized ground state and lower the energy further. In princi- 
ple, two type of excitations are available in a pyrochlore 
chain, either a magnon-like excitation produced by flip- 
ping a trimer state into its excited state or a pair of de- 
confined objects composed of a dimer plus a monomer. 
For our system with doubly degenerate ground state, the 
spinon-like excitations have lower energy. 

The deconfined excitations behave like domain-wall 
solitons which connects two spontaneously trimerized 
ground states. Explicitly, we represent an excited state 
with a dimer at site 2m- 1 and a monomer at site 2n is 
represented as 

5* (m, n) 

= ■ ■ - T (5to_i_i, Srn-1,2, Sm-l,^) m{Sm,l) 

T {Sm,2, Sm,3, <S'jn-|-l,l) • • • T {Sn-1,2, Sn-1,3, Sn,l) 
d {Sn,2, Sn,3) T {Sn+lS, S'n+1,2, <5'n+l,3) ' ' ' , 



where d{Si,Sj) = {\ai/3j) — \Piaj)) with a ^ (3 rep- 
resents a dimer. The corresponding momentimr-space 
wavefunction is 

l<m.<n<M 

The excitation spectrum can be calculated directly by 
using the above variational wavefunction. Because there 
exists no intrinsic mechanics responsible for binding the 
dimer and monomer together to form a bound state 
in a spontaneously trimerized system, it is reasonable 
to assume that the dimer and monomer are well sep- 
arated and they could be treated separately. Similar 
schemes have been used to evaluate the excitation spec- 
trum in the spin-sawtooth systeroi^i^. Under such an 
approximation, the excitation spectrum can be repre- 
sented as a sum of monomer part and dimer part, i.e., 
uj (fc,„, kd) ~ uj (km) + uj (kd) ■ The state ^' (n) is not or- 
thogonal with the inner product given by 



(*K)II*H) = 




thus 'i' (kd) has a nontrivial norm {'i' {kd)\\'^ {kd}) = 
4/(5 — 3 cos kd). With a similar scheme as solving the 
spectrum of excitations of the spin- 1/2 mode l^^i^^'^° , we 
get the spectrum for the dimer 

/, N 5 , 3 cos fed . 
Lo{kd) = -Ad (20) 

with Ad = 2 J. In the large J± limit, a monomer can only 
hop around in the left and right phases along the horizon- 
tal direction without breaking additional trimer singlet, 
therefore its excitation spectrum form a flat band, i.e., 
^ {km) — 0. When J± is close to J, the monomer actually 
can move around several corners, i.e., the monomer can 
be in the site of (m, 2) and (to, 3) either, therefore the 
wavefunction for a monomer has resonating structure at 
the mth triangle in the pyrochlore. The process of hop- 
ping from site (to, 1) to (to, 2) or (to, 3) accompanies with 
an energy cost of J± — J, and thus the excitation spec- 
trum is still a flat band for the homogeneous pyrochlore 
with J± — J. For the inhomogencous case, we take the 
monomer excitation as a three-site cluster consisted of 
three single monomers at sites (to, 1), (m, 2) and (m, 3), 
i.e., 

*(m) = -^[- ■ ■m{Sm,i) \ rn{Sm,2)--- 

H m{S.m,3) ■ ■■]■ 

Similarly, after considerable algebra, we get the spectrum 
for the monomer 

2 6 cos k 
ijj{km)^7;Am — (21) 

3 15 

with A„i ~ J± — J . 
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V. CONCLUSIONS 

We have proposed and studied a class of frustrated lat- 
tice which can be viewed as a ID strip of the pyrochlore 
lattice or a tetrahedron chain. For the general Heisenberg 
exchange model, we give an exact proof for the existence 
of the ground state consisted of the rung-dimerized state 
and the ground state of the decoupled chain. The phase 
diagrams of the spin- 1/2 and spin-1 tetrahedron chain 
are given and the phase boundaries are precisely deter- 
mined for the small-size systems. For both the spin-1/2 
and spin-1 systems, there; cixist thrcx; phascis. say, the fully 
dimerized singlets on the rungs plus a decoupled chain, 
a mixed phase with alternating spin singlet and state 
with total rung spin 2S on the rungs and a ferrimagnetic 
phase with long-range order, as the strength of vertical 
exchanges varies from infinity to minus infinity. We also 
studied the SU (3) spin-exchange model on the ID tetra- 
hedron chain, for which the ground sate is a double de- 



generate trimcrizcd state and the elementary excitations 
are fractionalized topological excitations. Our results in- 
dicate that the properties of the ground state for the 
pyrochlore systems with half-integer and integer spins or 
the systems with the same spins but different internal 
symmetries {SU{2) and SU{^) for spin-1 systems) are 
quite different. 
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